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ABSTRACT 

We report the results from a detailed analysis of an archival XMM- Newton ob- 
servation of the X-ray source XGPS-I J183251-100106, which has been suggested 
as a promising magnetic cataclysmic variable candidate based on its optical prop- 
erties. A single periodic signal of ~ 1.5 hrs is detected from all EPIC cameras on 
board XMM-Newton. The phase-averaged X-ray spectrum can be well-modeled 
with a thermal bremsstrahlung of a temperature kT ~ 50 keV. Both X-ray spec- 
tral and temporal behavior of this system suggest it as a eclipsing cataclysmic 
variable of AM Herculis (or polar) type. 

Subject headings: binaries: close — cataclysmic variables — stars: individual 
(XGPS-I J183251-100106, 2XMM J183251.4-100106) — X-rays: stars 



1. Introduction 

XMM-Newton Galactic plane survey (XGPS) has revealed a large population of X-ray 
sources (Hands et al. 2004). In order to identify the nature of the brightest sources detected 
in XGPS, an optical campaign has recently been carried out (Motch et al. 2010). In this 
process, three sources were identified as promising cataclysmic variables (CVs). 

The X-ray source designated as XGPS-9 (=XGPS-I J183251-100106) in Motch et al. 
(2010) is one of the newly found CV candidates in XGPS. Two optical objects are found 
within its X-ray error circle, while the fainter one (V ~ 23.3) is suggested to be the promising 
counterpart. Its CV nature was identified through optical spectroscopy which unambiguously 
shows the strong H, He I and He II emission along with a blue continuum (Motch et al. 2010). 
These features are known to be typical for a magnetic CV. 
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Apart from the optical identification, Motch et al. (2010) have also reported a brief X- 
ray analysis of archival XMM-Newton data obtained by the observations on 15 March 2002 
(Obs. IDs 0135741601, 0135744401). The authors found that an absorbed single component 
thermal plasma model is able to describe the data. However, as XGPS-9 located at an off- 
axis angle > 9' in these short exposures (~ 8 ks and ~ 5 ks respectively), the poor statistic 
and the degraded angular resolution leave the spectral parameters unconstrained. While the 
X-ray light curve of XGPS-9 (Fig. 8 in Motch et al. 2010) clearly shows the variability of 
this source, these short observations preclude the search for any periodic behavior of this 
system. 

In order to tightly constrain its X-ray properties, we have performed a follow-up inves- 
tigation of XGPS-9 with a deep XMM-Newton observation. The results of this analysis are 
presented in this paper. 

2. Observation & Data Analysis 

XGPS-9 was observed by the European Photon Imaging Camera (EPIC) on board 
XMM-Newton, which consists of two Metal Oxide Semiconductor CCD detectors (MOS1 
and MOS2) and a PN CCD detector, on 8-9 October 2009 (Obs. ID: 0605480101). While 
the MOS1/2 cameras were operated in Full Frame mode with a a temporal resolution of 
2.6 s, the PN camera was operated in Extended Full Frame mode with a temporal resolution 
of 199.1 ms. Medium filters were used to block optical stray light in all EPIC cameras. 

The aimpoint of this observation is RA=18 h 32 m 57.01 s Dec=-10°05'41.0" (J2000). Using 
tasks emproc and epproc of the XMM Science Analysis Software (XMMSAS version 11.0.0), 
we have reprocessed all the EPIC data with the updated instrumental calibration. We 
subsequently selected only those events for which the pattern was between — 12 for both 
MOS cameras and — 4 for the PN camera in 0.3 — 12 keV. We further cleaned the data 
by accepting only the good times when sky background was low for the whole camera (< 
2.4 counts/s, < 2.8 counts/s and < 5.7 counts/s for MOS1, MOS2 and PN respectively). 
After removing all events potentially contaminated by bad pixels, the effective exposures are 
found to be 55.9 ks, 56.6 ks and 41.3 ks for MOS1, MOS2 and PN respectively. 

We have also utilized the Optical Monitor (OM) data for our investigation. XGPS-9 
has been observed by OM in standard imaging mode with three different filters: U (effective 
wavelength 3440A), UVW1 (2910A) and UVM2 (2310A). The exposures of OM with U, 
UVW1 and UVM2 are 4.4 ks, 4.4 ks and 2.5 ks respectively. For data reprocessing, flat- 
fielding and source detection, we have utilized the metatask omichain. 
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This observation was originally intended for spectro-imaging investigation of the super- 
nova remnant Kes 69 (see Figured]), which has its X-ray properties poorly constrained (cf. 
Yusef-Zadeh et al. 2003, Bocchino et al. 2012). A detailed analysis of the remnant emission 
from Kes 69 will be published elsewhere (Seo et al. in preparation) as this is out of scope 
in this paper. In this observation, XGPS-9 is an serendipitous source located at an off-axis 
angle of ~ 5' (Figure [1]). With the aid of the XMMSAS task edetect-chain, we have run 
source detection individually on each EPIC data set. The mean X-ray position of XGPS-9 
determined from this observation is RA=18 h 32 m 51.51 s Dec=-10°01'05.0" (J2000) with a 
resultant uncertainty of 0.44" by combining the statistical errors inferred from each camera 
in quadrature. This X-ray position is marked as a black cross in Figure [H 

For timing analysis, the point source was extracted within a circular region of a 15" 
radius centered at the X-ray position from each camera, which corresponds to an encircled 
energy fraction of ~ 70%. 758 counts, 682 counts and 1277 counts were obtained from MOS1, 
MOS2 and PN respectively. All the photon arrival times are subsequently corrected to the 
solar system barycenter with the XMMSAS tool barycen by adopting the updated ephemeris 
JPL DE40£@ and the aforementioned mean X-ray position. 

In order to search for any periodicity, we have applied the techniques of both epoch- 
folding and the Lomb-Scargle periodogram (Scargle 1982) on the event list and the binned 
light curve respectively. The epoch-folding method was performed with the efsearch sub- 
package of HEAsoft. The peak value of the \ 2 spectrum with 50 bins for the arrival time 
of merged MOS and PN events is at P = 5337.7 ± 21.7s. The uncertainty of epoch-folding 
is estimated according to the empirical formula of Leahy (1987). On the other hand, we 
have binned all of the events into a light curve with 150 s resolution to calculate the Lomb- 
Scargle periodogram. The periodogram demonstrates a periodicity at P = 5353.0 ± 13.5s 
which is consistent with the result of epoch-folding. The uncertainty was estimated by 10 4 
times of Monte Carlo simulation. The results of the periodicity search with epoch-folding 
and Lomb-Scargle method are shown in Fig. [2j 

In order to identify the nature of this periodicity, we have also folded the photon events 
based on the period we detected. The 5337.7s period was adopted to fold the light curve as 
it is directly resulted from epoch-folding. In addition, we have divided the photons into soft 
band (0.3-2.5 keV) and hard band (2.5-12 keV) for computing the hardness ratio: (hard- 
soft) /(hard+soft). The folded light curves of both bands and the hardness ratio are shown 



4 We note that in comparing with DE405, the previous generation of planetary ephemeris, DE200, has a 
~ 300 m error in Earth's position. This corresponds to a timing error of ~ 1 /is which is negligible in our 
investigation of orbital period. 



-4- 



in Fig. [31 where the epoch zero is set at MJD 55113.0383. 

The minimum at phase zero is particularly interesting and leads us to further investigate 
it. This feauture resembles those resulted from the eclipses of X-ray emitting regions in polars 
(e.g. Mukai et al. 2003; Vogel et al. 2008). We notice there is residual X-ray emission at 
the minimum. To investigate if the residual is caused by the background, we estimated the 
count rates within a number of source-free circular regions of 15" radius around XGPS-9. The 
average background count rate estimated from the merged dataset is (3 ±2) x 10~ 3 counts/s 
which is consistent with the minimum of the folded light curve in 0.3-12 keV. This suggests 
the residual X-rays are likely the background contamination. Therefore, the feature can 
possibly be a total eclipse of the X-ray emitting region by the companion star. 

To further characterize this feature, we attempted to estimate the eclipse width and 
to obtain an ephemeris. In view of the relatively low statistics, these properties cannot be 
properly constrained by fitting the eclipse profile. Instead, visual inspection was adopted for 
the measurements reported in this paper. The folded light curves shown in Figure [3] suggests 
the ingress and egress are rather sharp which allow us to estimate the eclipse duration with 
the flat bottom minimum. This leads to an approximate width of 320 ± 53 s, where the 
uncertainty is estimated from half of the bin size in Figure [3j 

For defining the mid-eclipse phase, we adopted the mid-point of this feature in Figure [3j 
The timing of the X-ray eclipse can be described by the following linear ephemeris: 



T^id-edipsc) = MJD(TDB) (55113.0383 ± 6.2 x 10~ 4 ) + — — — — x N. (1) 

o04UU 

We have also checked these results by adopting different bin size (32, 64, 128 bin/period), 
we found that both the eclipse width and the ephemeris are consistent with the aforemen- 
tioned results within the tolerence of uncertainties. 

To investigate if the periodic signal comes from the synchronous rotation, we further 
examined the existence of any marginal periodic signal in a larger frequency range using 
the x 2- test with 32 bins and Fourier resolution. Because the spin period of a typical white 
dwarf is about hundreds of seconds (cf. Ritter & Kolb 2003) and the signal of orbital 
modulation seriously contaminates the periodogram when the searching range is larger than 

2.5 xlO -3 s _1 , we only considered a periodic signal of spin within 100-400 s. Within this 
selected range, the largest peak value can be found at 164.7(1) s with the chance probability of 

9.6 x 10~ 5 . However, this result is not significant if we took into account the number of trials 
in the searching range. We also tried the decomposition methods (e.g. the empirical mode 
decomposition proposed by Huang et al. 1998) to resolve another strong signal. However, 
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we did not yield any other possible periodicity except for the binary modulation. Therefore, 
we conclude that no spin period can be detected from the current XMM-Newton observation 
and it may give a hint of the synchronized spin and orbital periods. 

Motch et al. (2010) have briefly reported the X-ray spectral properties of XGPS-9. 
They found that a one-component thermal plasma spectrum with a column absorption of 
Nh ~ 8 x 10 21 cm -2 can provide a reasonable fit to the data. However, due to the limited 
photon statistic, the plasma temperature and the X-ray flux in their study remains to be 
unconstrained. With a much longer integration time of the data adopted in our investigation, 
we are able to provide a tighter constraint on its spectral properties. For the spectral analysis, 
we have extracted the spectra within a circular region with a 30" radius in each MOS camera 
and a 25" radius in PN camera around the X-ray position respectively. This corresponds to an 
encircled energy fraction of ~ 85% in all EPIC cameras. Spectra extracted from each camera 
were grouped so as to have at least 50 counts per spectral bin. The background spectra were 
sampled from the low count circular region of 30" radius centered at RA=18 h 32 m 55.72 s , 
Dec=-10°02 05.87 (J2000) in each camera. Response files were computed by using the 
XMMSAS tasks rmfgen and arfgen. For the spectral analysis, we used XSPEC (version 
12.6.0) with x 2 statistics adopted for all the fittings. All the quoted uncertainties of the 
spectral parameters are la for 1 parameter of interest. 

Based on the optical properties, Motch et al. (2010) suggested XGPS-9 a magnetic 
CV. To explore this scenario, we examine its X-ray spectrum by fitting an absorbed thermal 
bremsstrahlung model, which accounts for the X-ray emission from the shock-heated gas in 
the accretion column. We found that this single component model can describe the data 
reasonably well with a goodness-of-fit of x 2 =49.40 for 54 d.o.f.. The spectral fit yields a 
column density of = 7.6l '7 x 10 21 cm" 2 , a plasma temperature of kT = 46 ± 10 keV and 
an unabsorbed flux of f x = 3.3^02 x 10~ 13 erg cm" 2 s" 1 in the energy range of 0.3 — 10 keV. 
The comparison between the best-fit model and the data is shown in Figure HI We have also 
computed the confidence contours to investigate the relative parameter dependence between 
the plasma temperature and the column absorption, which are plotted in Figure [5j 

As X-ray spectrum of CVs can possibly show the presence of line emission, we have also 
examined the observed spectrum of XGPS-9 with XSPEC model MEKAL which is a code 
that models the plasma in collisional ionization equilibrium with the line emission built in. 
We found that the fit yields a comparable goodness-of-fit with the thermal bremsstrahlung 
model (x 2 =50.11 for 54 d.o.f.). The best-fit spectral parameters are also similar (A^h = 
7.5 ± 0.7 x 10 21 cm" 2 , kT = keV). 

For magnetic CVs, part of the thermal bremsstrahlung emission can possibly be Comp- 
ton reflected by the white dwarf surface (see Matt et al. 1998). This so-called Compton 
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reflection component is expected to be accompanied with a fluorescent Ka iron line (de 
Martino et al.2008). For searching evidence for the line feature, we have added a Gaussian 
component to the thermal bremsstrahlung model and fixed the line energy at 6.4 keV. To 
begin with, we have allowed both normalization and the line width as free parameters. How- 
ever, we found that the line width cannot be constrained (i.e. the fitted width is essentially 
zero). This led us to fix it at a reasonable value of a = 100 eV (see Balman 2011; de 
Martino et al. 2008). This results in a best-fit parameter set of N H = 7-8-o' 9 $ x 1q21 cm~ 2 , 
kT = 36 ± 10 keV, and a line flux of / 6 . 4 kc v = 3.0+^ x 10~ 7 photons cm~ 2 s" 1 with 
X 2 =49.11 for 53 d.o.f.. In view of the large uncertainty of fitted line flux and the insignif- 
icant improvement of the goodness-of-fit, we conclude that there is no compelling evidence 
of Fe Ka line can be found in this observation. We report a la upper limit on the line flux 
of feA kcv < 7.5 x 10 -7 photons cm~ 2 s _1 . 

While a thermal bremsstrahlung/MEKAL can model the spectrum well, as mentioned 
by Motch et al. (2010), the spectral hardening can possibly be resulted from the reflection 
component and/or warm absorption (see Schwarz et al. 2009; Staude et al. 2008). There- 
fore, the plasma temperature inferred from this simple model can be unphysically high. We 
have also investigated if any other single component model can provide any acceptable phe- 
nomenological description. We notice that a simple absorbed power-law can also result in 
a comparable goodness-of-fit (x 2 =49.83 for 54 d.o.f.), and hence we cannot distinguish it 
from the thermal bremsstrahlung/MEKAL model statistically. The power-law spectral fit 
yields A H = 7.9+JJ x 10 21 cm -2 , a photon index of T = 1.3 ± 0.1 and an unabsorbed flux 
of f x = 3.4+g;7 x 10~ 13 erg cm" 2 s -1 (0.3 - 10 keV). The spectral steepness and the flux 
are fully consistent with those inferred from a brief analysis limited to the hard band (i.e. 
3 — 10 keV) as reported recently (cf. Tab. 1 in Bocchino et al. 2012). 

We have checked the robustness of all the spectral results quoted in this paper by 
incorporating background spectra sampled from different source-free regions around XGPS- 
9. We found that the spectral parameters inferred from the fittings with different adopted 
background are consistent within la uncertainties. 

We have also attempted to constrain the possible contribution from the soft component 
by using OM data. Soft component from a CV can be originated from the optically-thick 
blackbody like emission of the white dwarf surface. The temperature can range from few 
tens to hundred of eV (Evans & Hellier 2007). Searching for the optical/UV counterpart 
in OM data results in non-detection in all three bands, which place a limiting flux density 
of < 2.3 ^Jy, < 2.7 /iJy and < 7.8 //Jy in U, UVW1 and UVM2 respectively. We adopted 
the column density inferred from the X-ray spectral fitting (i.e. 7.5 x 10 21 cm~ 2 ) to per- 
form the extinction-correction (cf. Predehl & Schmitt 1995; Cardelli et al. 1989). With 
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the de-reddened limiting flux densities, we place an upper bound for any soft component 
contribution by assuming a low temperature white dwarf with kT = 10 eV (see Figure |6]). 
The blackbody gives a limiting flux of < 5.4 x 10~ 18 erg cm -2 s _1 in 0.3 — 10 keV, which is 
< 0.002% of the total X-ray flux in this energy range. 

3. Summary & Discussion 

Utilizing the archival XMM-Newton data, we have detected the periodicity from a new 
CV candidate XGPS-9. Through a detailed timing analysis, only a single periodic signal of 
~ 1.5 hrs can be found. This suggests a possible scenario of synchronous rotation, which 
favors the interpretation that the system is a magnetic CV of polar type. 

For a polar, the white dwarf is accreting matter from its late-type main sequence com- 
panion via Roche-lobe overflow which will be driven by the strong magnetic field of the 
white dwarf (with a surface magnetic field strength of few tens of million Gauss) to form a 
quasi-radial flow towards the magnetic pole. The accretion flow will eventually form a strong 
shock above the stellar surface. The shock-heated plasma will be cooled via bremsstrahlung 
radiation in X-ray regime with typical temperature of few tens of keV, which is consistent 
with the best-fit plasma temperature (i.e. kT = 46±10 keV) inferred in our spectral analysis. 

It is interesting to compare the temporal behavior of XGPS-9 and 2XMMp J131223. 4+173659, 
which is a eclipsing polar serendipitously identified in 2XMM catalog (Vogel et al. 2008). 
The properties of these two systems are remarkably similar. First, both objects show a nar- 
row and deep minimum in their phase-folded light curves (i.e. depression in phase 0.97 — 1.03 
for the case of XGPS-9). This feature is interpreted as the eclipse of X-ray emission in the 
polar region of the white dwarf when the companion star passes through our line-of-sight 
towards this region. Also, the broad depression (i.e. phase ~ 0.2 — 0.6 in Fig. EJ), which 
is likely due to the self-occultation by the white dwarf itself, is observed in both system. 
The coexistence of these two features indicates that the spin of the white dwarf and the 
binary orbital periods are synchronized. Furthermore, a shallow minimum in the soft band 
light curve (in phase ~ 0.7 — 0.9 in Fig. [3]) is found before the eclipse. When the X-ray 
emitting region is obscured by the accretion stream, a pre-eclipse dip can be resulted from 
the increased photoelectric absorption which makes the dip more prominent in the soft band 
(see Fig. [3J. 

We further attempted to place a constraint on the orbital inclination with our estimated 
eclipse width (i.e. ~ 320 s). Assuming the companion star fills its Roche lobe and X-ray 
emitting region is close to the surface of the white dwarf, the duration of the total eclipse of 
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the white dwarf can be expressed as a function of mass ratio and the inclination (cf. Fig. 2 
in Home 1985). The eclipse width corresponds to a phase interval of A(j) ~ 0.06 which we 
considered as an upper limit of the eclipse width as a longer observation might reveal the 
falling/rising edges in the eclipse profile. We further assume that the companion star is a 
M-dwarf with mass M* > 0.1 M . For a typical white dwarf (i.e. M w d ~ 0.6 M ), this 
implies a mass ratio of q = -ffi- > 0.17. With these assumptions, the upper bound of the 
orbital inclination can be placed at < 80°. 

For further investigation of XGPS-9, dedicated optical/infrared observations will be 
useful for confirming its source nature. A differential photometric study can enable a com- 
parison with the X-ray light curve reported in this paper. This can provide information 
for further constraining the accretion geometry. Also, the free electrons in the shock region 
can spiral in the magnetized ionized gas. This can lead to the emission of cyclotron radi- 
ation in infrared regime. Therefore, infrared spectroscopy and polarimetry are encouraged 
for determining the magnetic field of the white dwarf directly. We would also like to point 
out that the current data cannot provide a conclusive constraint on the plasma temperature 
of a thermal bremsstrahlung spectrum at a few tens of keV, which is far above the energy 
coverage of XMM-Newton. This is reflected by its large relative uncertainty reported in 
this work (i.e. > 20%). For a more constraining spectral results, hard X-ray observation 
is certainly needed. As XGPS-9 resides in rather complex region, the recently commenced 
mission NuSTAR and the upcoming one (e.g. Astro-H), which equipped with hard X-ray 
focusing optics, will provide the desirable instruments for the further investigations of this 
interesting eclipsing polar. 
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Fig. 1. — X-ray image of a 18' x 18 FOV toward the supernova remnant Kes 69 with the 
MOS1/2 and PN data merged. XGPS-I J183251-100106 is the bright serendipitous source 
marked by a black cross. The dashed circle represents the background region adopted for 
the spectral analysis (see text). 
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Fig. 2. — The result of epoch-folding period search (upper) and the Lomb-Scargle peri- 
odogram (lower). The dashed line in the Lomb-Scargle periodogram is the 3a significance 
level. 
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Fig. 3. — The folded light curve of soft band (upper), hard band (middle), and the corre- 
sponding hardness ratio (bottom). The unit of y— axis for the first two panels is count s _1 . 
We have adopted a resolution of 50 bins per cycle. 
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Fig. 4.— Energy spectrum of XGPS-I J183251-100106 as observed with the PN (upper 
spectrum) and MOS1/2 detectors (lower spectra) and simultaneously fitted to an absorbed 
thermal bremsstrahlung model {upper panel) and contribution to the x 2 fit statistic (lower 
panel). 
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Fig. 5. — la, 2a and 3a confidence contours for the thermal bremsstrahlung model fitted to 
the X-ray spectrum of XGPS-I J18325 1-100106. The x 2 values for the minimum and each 
contour are given at the top of the figure. 
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Fig. 6. — The dash line represents the upper-bound of possible blackbody contribution at 
kT = 10 eV as constrained by the extinction-corrected limiting flux densities deduced from 
the OM data (triangular symbols). 



